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ABSTRACT 

We present a new method for determining the column density of erupting filament material using state-of-the- 
art multi- wavelength imaging data. Much of the prior work on filament/prominence structure can be divided 
between studies that use a polychromatic approach with targeted campaign observations, and those that use 
synoptic observations, frequently in only one or two wavelengths. The superior time resolution, sensitivity and 
near-synchronicity of data from the Solar Dynamics Observatory's Advanced Imaging Assembly allow us to 
combine these two techniques using photo-ionisation continuum opacity to determine the spatial distribution 
of hydrogen in filament material. We apply the combined techniques to SDO/AIA observations of a filament 
which erupted during the spectacular coronal mass ejection on 2011 June 07. The resulting "polychromatic 
opacity imaging" method offers a powerful way to track partially ionised gas as it erupts through the solar 
atmosphere on a regular basis, without the need for co-ordinated observations, thereby readily offering regular, 
realistic mass-distribution estimates for models of these erupting structures. 

Subject headings: methods: data analysis radiative transfer Sun: activity Sun: coronal mass ejections (CMEs) 
Sun: filaments, prominences Sun: UV radiation 



1. INTRODUCTION 

The mass and structure of prominences/filaments are of in- 
terest in several aspects of the study of the corona. Since the 



launch of Hinode and its Solar Optic al Telescope (SOT; |Ko- 
|sugi et al.|2007||Tsuneta et al.|20 08), observations of the dy- 



namic structures m prominences with SOT data (jBerger et al. 



|2011| ) have underlined the difficulties of explaining promi 
nenc e dynamics. Development of radiative transfer models 



(e.g. Labrosse & Gouttebroze|2004[ |Anzer & Heinzel 2005 ) 



requires an understanding of the fine- scale mass structure 
and radiation field, while oscillations in prominences can be 
used for diagnosis of the magneti c field streng th only if their 
mass can be accurately estimated ( Oliver 2009 , and references 
therein). Models th at seek to replicate the initial conditions of 
filament eruptions ( [Low et a l. 2003 ), and th ose that address 
the ensuing propagation into the heliosphere dDeForest et al. 



ag propagation into me heliosphi 
1 2012]) must take into account t he effects of gravi ty ( |Spicer 



Fani2006l and composition ( |Kilper et al.| r2009)' on these 



concentrations of material that are at least an order of magni- 
tude denser than the surrounding corona (e.g., Bommier et al. 
[T994HJejcic & Heinzel|2009l ). 

Tne most established methods of measuring mass distribu- 
tions use Thomso n scattering of photosp heric light by free 
coronal electrons ( Hundhau sen et al.|1994| ) to estimate the to- 
tal amount of mass in promin ences and (interplanet ary) coro- 
nal mass ejections ((I)CMEs; Vourlidas & Howard 2006, and 
references therein). More recently, investigators have used 
the amount of material evacuated in the dimming of the EUV 
corona to estimate the amount of mass lost to a CME/ICME 



(Harr ison & Lyons |2000[ [Aschwand en et al.|20 09), showing 
consistency with parallel observations from Thomson scatter- 
ing estimates. 

At wavelengths shorter than the Lyman series limit (912 A), 
cool plasma is optically thick due to photo-ionisation, princi- 
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pally of neutral hydrogen and neutral or singly-ionised he- 
lium. The product of abundance and cross-section for photo- 
ionisation of other species is negligible compared with these 
three. This continuum absorption is used in studies of the 
inters tellar medium (see |Vennes et al.||1993[ |Fruscione et al.| 
1994, and references therein) to calculate the column density 
(N( j) of these species along the line of sight, where i is the 
element, and / the ionisation stage (/ G • • ■ ,Z}), since 
the total opacity at a given wavelength, t(A) is given by 



(i) 



where Ou is the cross-section for absorption due to photo- 
ionisation for species U. This allows the column mass M 
along the line of sight to be estimated as: 



(2) 



The effect was also note d in the solar transition region by 
|Schmahl & Orrall| ( |1979| ) as evidence of H I and He I - and 
therefore partial ionisation - in the chromosphere and transi 



tion region. (See also the work of |Orrall & Schmahl 1976 
|Kanno|1979l|Kanno & Suematsu|1982f ~ 



The most frequent investigations of column-density struc- 
ture on the Sun are those that estimate how much light is 
attenuated by a filament, - i.e., its opacity - and, so, total 
mass along the line of sight in each spatial resolution ele- 
ment (Equation [2]). Since both band-pass imager and slit- 
spectrometer data of the Sun exist in the EUV, these inves- 
tigations have consisted of various combi nations of images a t 
a single wavelength or wavelength band ( Gilbert et al. 2006 ), 
images at multiple wavelengths ( Golub et al. 1999; Engvold 



|et al. |200"l} , spect ral line m easurements in the EU V and (o f- 
te n) Ha data (e g , |Schmieder et al.| 19991 |Mein et al.|2001b . 

Gilb ert et al.| ( |2005| |2006| ) have developed and applied a 
technique for estimating the mass of limb-crossing promi- 
nences using EUV images at 195 A only. This so-called 
"spatial-interpolative" approach estimates the background 
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emission and prominence opacity. An interesting point is that 
they make use of the fact that erupting prominences eventu- 
ally reveal something close to the original background emis- 
sion, albeit at a later point in time. This allows a "temporal- 
interpolative" approach, to which we return later. Though not 
strictly "monochromatic", we use the term in this article to 
refer to all single- wavelength/-passband methods. 

A logical extension of the above method is to use the in- 
formation from several EUV wavelengths, over a sufficiently 
broad range, so that the cross-section for photo-absorption 
varies while other variables are assumed to hold constant. If 
all the associated radiation escapes from the corona, where 
the pressure scale-heights are large, then this approximation 
can be made. Such "polychromatic" approaches include: 



2010) where, for coronal lines, there is lower emissivity in 
both the cavity and the filament itself. This effect is referred 



those which are (coronal) image-based (Golub et al. 1999 



Gilbe rt et al.|2011|); those that use coronal li ne spectroscopy 
(Kuc era et al.||1998| |Anzer & Heinzel||2005| ); and those that 
include emi ssion from the t ransition regi on and/or ch romo- 
sphere (e.g., Schmieder et al. ||1999} |Mein et"aL]|2001 |), often 



interpreted with the aid of non-LTE radiative transfer mod- 
elling to gain further insight into the detailed ionisation and 
mass structure of a fila ment. Motivated to con strai n filament 
mode l s, for example, | Anzer & H einzel (2003) and |Schwartz 



|et al,| P004| ) used a radiative model (Heinzel et al.||2003| rto 



infer the height of filament material above the solar surface, 
and to differentiate between morphological filament models. 



The third o f the listed techniques is still developing. [Labrosse 
et al.| ( |20l"T] ) most recently a pplied it to Hinode EUV Imaging 



Spectrometer spectra (EIS; Culhane et al. 2007), while Gilbert 
|et al. | d20 1 1|> exp loited SoHO imaging and spectroscopic data 
(Domingo et al.||1995|). We ref er the reader to the recent re- 
view by Labr osse et al.| ( |2010| ) for an excellent treatment of 
these and other prominence diagnostics. 

A distinct disadvantage of the spectroscopic approaches, 
pragmatically speaking, is that contemporary EUV spectrom- 
eters have a limited field of view, so that the observations must 
be targetted. Ha measurements, on the other hand, are typi- 
cally ground-based and so are susceptible to weather-imposed 
limitations. As a result, these methods have not been used on 
a regular basis. 

A persistent unknown in much of this work is the back- 
ground radiation field, here denoted 4, that is attenuated by 
the filament material. The background corona is clearly vari- 
able on short and long spatial and temporal scales, but a view 
of it is blocked by the material under investigation. This 
has been treated by interpolating bet ween points in the back- 
ground on either side of cool masses ( Kuce ra et al.|19 98), with 
additional consider ation of the profileoT Background corona l 
emission off-limb (Gil bert et alJ|20TT| |Golub et al.|[T999] ). 
However, because of the highly structured nature of the so- 
lar atmosphere, a better background model can be obtained 
by observing coronal emission in the absence of the absorb- 
ing structure, such a s before and/or after a filament eruption. 
Gilb ert et"aLl ( 2005| ) used both a "temporal-interpolative" ap- 
proach and a "spatial-interpolative" approach to measure the 
column depth of a prominence which erupted from the south- 
east solar limb; their method was later applied to a larger sam 



le of erupting and non-erupting prominen ces ( [Gilbert et al. 
2006| ). Similarly, |Kucera & Landi| ( |2008| ) use a temporal- 
mterpolative approach in analysis of data from an erupting 
prominen ce structure taken with the SOHO/SUMER spec- 
trometer dWilhelm et al.|1995| ). 

It is worth pointing out that photometric measu rements are 
complicated by the presence of filament cavities ( Gibson et al. 



to as "vol ume blocking" or "emissivity blocking" (Heinzel 



|et al. |2003| ), and simply refers to a lack of emission measure. 
The amount of attenuated radiation emerging from behind the 
filament is difficult to disentangle from radiation emitted in 
the foreground, but a promising alternative to the interpola- 
tive appr oaches is to use th e X-ray transparency of filament 
material ( Anzer et al. 2007) and assume that the EUV back- 
ground emis sion can be scaled to match the spatial soft X-ray 
distribution ( [Heinzel et al.||2008| ). Whilst this, too, is an as- 
sumption, it does address the issue of reduced emission in the 
cavity. 

Since the launch of NASA's Solar Dynamics Observatory 
(?) incorporating the Advanced Imaging Assembly (AIA; ?), 
we have access to much higher-cadence observations of the 
solar corona, with data simultaneous or near- simultaneous in 
several EUV wavelengths. This allows us to apply a temporal- 
interpolative approach to measure the total background and 
foreground emission around cool-material structures in appar- 
ent motion across the Sun, often at several points close in time 
to the transient absorption taking place. 

In this article, we describe the application of two methods 
- one monochromatic, the other polychromatic - to AIA ob- 
servations of cool material returning from a spectacular fila- 
ment eruption. Although the eruption on 2011 June 7 is suc- 
cessful, producing a CME, some filament material returns to 
the solar atmosphere and surface, and our aim is to infer the 
density of material from its obscuration of the background 
corona. We describe two approaches to achieving this aim. 
The first approach calculates the opacity of the material from 
data taken in a single AIA fi lter, an approach also used in [van] 
|Driel-Gesztelyi et aT] ( [20 1 2| ) but described here in detail. The 



second approach uses arguments first made by Ku cera et al. 
(1998) to separate the column density from the product of fill- 
ing factor and geometrical depth using observations made at 
multiple wavelengths. In both cases, we use properties of the 
opacity due to H I, He I and He II photo-ionisation continua to 
give a more powerful lower-limit estimate of the total hydro- 
gen column density along the line of sight, using the vastly 
improved cadence and sensitivity of AIA data. Lastly, we 
compare the results of these approaches. 

2. OBSERVATIONS 

The observations used in this work were recorded by AIA 
on 2011 June 7, in the period 06:30 - 08:00 UT. The region 
of interest is around NOAA active region 11226, which pro- 
duced an M2.7-class X-ray flare and a spectacular filament 
eruption, in which most of the filament mass appears to have 
returned to the Sun. The returning filament mass is the sub- 
ject of this study, where we estimate the column and volumet- 
ric densities. We focus on two primary targets in which to 
diagnose the total hydrogen column density, A/r. 

Target 1 is a large concentration of returning material that 
appears to fall unhindered by solar magnetic field, suggesting 
high plasma j3, consistent with a low ionisation degree. The 
material is shown in Figure [TJ). 

Target 2 is an area to the east of the active region, where re- 
turning filament mass arrives at what resembles a Y-shaped 
coronal magnetic nu ll-point (Figure [6j see also |van D riel- 
|Gesztelyi et al.|2012 ). Material that appears dark in the EUV 
range enters and leaves this area and we estimate the density 
of one of the departing concentrations. 

In order to correctly estimate the measurement errors in our 
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data, we model the photon noise using the calibration curves 
distributed in the AIA branch of the SolarSoft IDL software 
library (SSW; |Freeland & Handy|1998| ). 

3. METHOD 

In this article, we assume that the dominant process which 
removes photons from the line of sight is photo-ionisation. 
We construct the cross-sections for photo-ionis ation, <T; / ( A ) , 
using the analytical approximations given by Verner et al. 
( 1996). We have also considered the effect of He I autoion- 
isatio n resonances, u sing Fano profile parameters given by 
Rum ph et al.| ( 1994] ). Although <7 He i can be enhanced by 
as much as a factor of 10 in these resonances, the effect is 
very narrow-band (full width at half-maximum < 0.5 A) when 
compared with the width of the AIA bandpasses (FWHM of 
AIA 193 channel = 6.3 A), so we disregard the effect of these 
resonances in calculating He I opacity. Additionally, we as- 
sume that there is no emission from the prominence material 
in the wavebands observed (we return to this issue briefly in 
Section [5}. 

As noted by Daw et aL] ( |1995 ), 0" H e i and <7 H e II are simi- 
lar in value (for A < 227 A). Here, however, we draw atten- 
tion to the fact that the cross sections of the first three species 
are very similar when weighted by elemental abundance, i.e., 
AiOu (Figure [2]). We follow the conv ention that Ah is unity , 
and take the value An e = 0.0851 from |Grevesse et aL] ( |2007| ). 
This lets us make the approximation 



*He = ^He I ° He I + #He II a He II 
-(A^Hel+^He II) ^ He II 
~A He AfH0He II 



(3) 



and since T H i(A) « TH e (A), for A < 227 A, x tot w 2Tn e . Al- 
lowing for the fact that some helium may be fully ionised, we 
can then estimate the total hydrogen column density along the 
line of sight as 

Nr > - . t>0 ' ■ (4) 
ZA He 0He ii 

The deviation of the abundance-weighted cross-sections of 
helium from that of H I is, in fact, less than a factor of 2 (see 
shaded box in Figure^) for all the wavelengths of AIA data 
analysed here. Because the above approximation is expressed 
in terms of the helium cross- sections, it allows us to more ac- 
curately gauge the total hydrogen column density for weak 
degrees of He II ionisation, independent of the ionisation de- 
gree of hydrogen, which fully ionises at temperatures much 
lower than does helium. Accordingly, and considering the in- 
equality sign in Equation [4] our estimate of A/r will be a lower 
limit; the lower the ionisation degree of helium, the more this 
value is likely to represent an accurate value rather than just a 
lower limit. 

3.1. Polychromatic measurements 



Kucera et al . | ( j 1 99 8) write the following expressions for the 
observed EUV intensity at a given point on a filament, con- 
sidering the effects of a pixel filling factor along the line of 
sight, /: 



obs ' 



where 



-I b [fe- T + (1 
-Io-fI b (l-e 



"/)] 
■ T ) 



/ 



(5) 



(6) 



This allows Equation[5]to be rewritten as: 

!_^ =/ £ (1 _ (7) 
*o 4) 

= G(l-e~ r ) (8) 

where the factor G combines the two unknown geometrical 
factors of the fraction of emission behind the absorbing mate- 
rial, 4 /To and the filling factor that de scribes how much of the 
pixel area is filled by this material. | Anzer & H einzel (2005 ) 
note that the former term can potentially be appreciably less 
than unity, particularly when the absorbing material is low in 
the atmosphere. The quantity on the left-hand side of Equa- 
tion[7]is an observable quantity, denoted the absorption depth, 
d. 



d=\ 



lobs 

d(X) = G[l-e 



-t(N b ; A)i 



F(G,N U ; X)=d(X) 



(9) 
(10) 

(11) 



In this case, d{X) can be fitted to a function F if there are in- 
dependent observations at a sufficient number of wavelengths. 

The above method is appropriate where we have simultane- 
ous observations of the same target. This also holds true for 
ft£tfr-simultaneous observations, provided that the target does 
not change shape or density appreciably between measure- 
ments. Thanks to the cadence of AIA, this is true in the case of 
Target 1. We therefore cross-cor relat e images of monochro- 
matic opacity, r(A) (see Section [3^2] ) to ensure accurate co- 
registration, since this property ought to scale largely with 
t(A) between measurements at different wavelengths. 

In this analysis, we use a Levenberg-Marquardt least- 
squares minimisation al gori thm to find the best fit to the mea- 
sured d(X) in Equation [10] An example fit is shown in Fig- 
ure [3] Since there are two free parameters of F, we require 
observations in at least three wavelengths to constrain them. 
We use data from the 94, 131, 171, 193, and 211 A channels. 
Although the images are near- synchronous rather than being 
truly cotemporal, the cadence of the observations is such that 
the form of Target 1 does not change appreciably between im- 
ages at these wavelengths. We disregard data from the 304 A 
channel as this emission is optically thick, while we are un- 
able to use data from the only other EUV channel (335 A) 
because of what appear to be stray- or scattered-light effects, 
exacerbated by the extremely bright flare emission from the 
active region where this eruption originates. The remaining 
AIA filters measure wavelengths that are longer than the Ly- 
man series limit. Figure [4] shows the resulting best-fit values 
of A^h and G for Target 1 . 

A difficulty in selecting the background area for this target 
is that for most of the duration of this event, many other con- 
centrations of erupted filament material cross the field of view 
in various directions, and only in the images taken around 
06:59:36 UT is there a clear view of the solar background 
emission - strictly, Iq - without attenuation. We therefore 
choose frames at (or nearest to) this time as our estimate of 
/o(A). 

3.2. Opacity measurements at a single wavelength 

The filament mass concentrations in Target 2 move, and de- 
form, much more rapidly than Target 1, so that it is not possi- 
ble to convincingly co-register images taken at different times 
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a) SDO 171 : 06:59:36UT 



b) SDO 171 : 07:03:36UT 
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Figure 1. a) Portion of an AIA 171 A filter image which is used as the model unattenuated image (7o) for Target 1, a concentration of filament material that has 
failed to escape during the CME under study, b) Target 1 itself, which has now fallen further towards the Sun. Dashed and solid intersecting lines refer to the 
positions of profiles taken through Nu shown in Figure|5] (See the electronic edition of the Journal for a color version of this figure, where it is also available as 
an MPEG animation.) 



(a) Abundance-weighted cross-sections for photoionisation 
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Figure 2. (a) cross-sections for photo-ionisation, AfGij, for the three prin- 
cipal contributors in the short EUV; (b) ratio of A/O/y to AhCT h j. The grey 
shaded box indicates (horizontally) the range of EUV wavelengths shorter 
than the He II ionisation edge that are observed by AIA, as well as the range 
within which the abundance-weighted cross-section can deviate by a factor 
of < 2 from that of H I (vertically). 



and wavelengths. This being the case, we cannot use near- 
simultaneous measurements of intensity reduction at different 
wavelengths to fit Equation [10| However, since we can no 
longer estimate G independently of Nh, we can simply set 
G = 1 and accept that we will make an underestimate of the 
column density in this way, since this implies a unity filling 
factor, and that all EUV emission Iq is behind the erupting fil- 
ament material. In this case, we would rewrite Equation 10 



for a single wavelength as simply 

T(A)>ln/ 0> A-ln/ ofa> A- 



(12) 



Using this monochromatic opacity, we are again able to cal- 
culate a map of the lower limit to this quantity, and of the re- 
sulting column density show in Figures H]& |6j: . Note that this 
type of opacity map is used co-register the nea r- simultaneous 



images at different wavelengths in Section 3.1 

The area covered by the field of view in Figure [6] is, fortu- 
nately, unobscured by returning filament material until around 
06:50 UT. There is a pronounced dimming around this region 
following the transit of a coronal wave, but the period between 
06:44 and 06:48 UT lies between the passage of the wave and 
the arrival of falling filament material. Therefore, we model 
the background emission, Iq, as the mean intensity of images 
taken during this interval (Figure [6^). Because we calculate 
Iq, in this case, from a time-average of these data, we estimate 
the measurement error on the background at each pixel as the 
quadrature sum of the photon noise from each image and the 
standard deviation of intensity in that pixel in time. 

4. ANALYSIS 
4.1. Polychromatic method 

We define the edge of Target 1 as the level where G = 0.6 
(Figure ffl. This level is set arbitrarily, but seems to match 
the visual edge of the target well in the portion of the FOV 
which was unobscured in Iq. The restriction of our analysis 
to the area indicated in Figure [4] reflects the fact that some 
portions of the image used as Iq were already obscured by ab- 
sorbing material. The animation that accompanies Figure [4] 
in the online edition of the journal shows that other, fainter 
filament material moves between the times indicated between 
Figures]?^ (constituting Iq) and |4j) (Jobs)* revealing an elon- 
gated bright structure behind Target 1. This would tend to 
cause an underestimate of Nu in the affected pixels, since the 
bright structure was not accounted for in Iq, so that the drop 
in intensity would in fact be slightly larger. 

The detailed map of best-fit column density (Figure |4^) 
shows a concentration of material towards the lowest part of 
the target, with column densities of Nu ~ 10 20 cm -2 . We 
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Figure 3. Example of the result of fitting observed d (X) at all AIA wave- 
lengths below 227 A. Data points are taken from a measurement inside the 
large mass shown in Figure (Th Large circles with error bars: d(X). Small 
circles joined by solid line: values of d due to the best fit to Equation [To| 
where the hatched envelope represents the bounds of the best fit ±la. (See 
the electronic edition of the Journal for a color version of this figure.) 

note that these values are rather larger than those previously 
reported in either erupting or quiescent f ilaments (jPenn 2000 ; 
Gilbe rt et al.|2 005 ; Labross e'et al.|2010] >. The numeric values 
and their associated 1<7 errors along two slices are shown in 
Figure [5] While the variance along a slice in is noticeable, 
the corresponding profiles of G are much less variable. In 
fact, the value of G seems remarkably stable across the target: 
given that Target 1 is likely to have reasonably uniform hjh 
on larger scales, its small-scale uniformity indicates that the 
pixel filling factor, / may be rather smooth. G also lies close 
to unity, which indicates that /, too, must be rather high. The 
absolute value of 4 /To is tricky to determine (without using 
stereoscopic data), but the high value of G indicates that the 
mass must be several emission scale heights into the corona. 

We should point out that the values along the slices include 
all the values fitted in the FOV shown in Figure |4j not only 
those within the target shown, which leads to the more wildly 
varying fitted values - from outside the G = 0.6 contour - 
being plotted at the edges of Figure [5^. 

In order to estimate mean volumetric (number) densities in 
this structure, we divide A/h by a characteristic path length 
through the material. We choose the FWHM of G (14") 
as our indicator, equivalent to a path length of / = 10 9 cm. 
The estimates of are therefore dependent on the value of 
/ being constant; although this is almost certainly an over- 
simplification for any given point in the target, it serves as a 
guide for estimating the total hydrogen density, which reaches 
nn > 10 10 cm" 3 . 

4.2. Monochromatic method 



As stated in § 3.2 the primary object to which we apply our 



monochromatic method is Target 2 (Figure [6p). In this case, 
we have used data from the AIA 193 A channel because of 
its s uper ior signal-to-noise ratio. In contrast to our approach 



m 



4.1 we do not isolate any one piece of material in the 



seen to be moving eastward (negative x direction) from the ap- 
parent Y-point. The positions of these profiles are indicated in 
Figure [6]:, and the corresponding values along each are shown 
in Figure [7] It is noticeable that Nu varies more smoothly 
across this target than in Target 1 (Figure [5]). Although the 
nature of the monochromatic analysis is a mapping from a 
more smoothly-varying intensity to T (rather than a fit of two 
free parameters t o mu ltiple data points), the background being 
time-averaged (§ 3.2) may also contribute to this smoothness. 

More generally, trie amount of material contained in these 
thinner 'threads' is still substantial, with lower limit values in 
the range A/h > 10 18 - 10 19 cm" 2 , similar to th ose seen in 
more extended prominences ( [Kucera et al.|1998| . In order to 
estimate the volumetric density tin, we again assume a repre- 
sentative path length equivalent to the FWHM of the structure 
in the plane-of-sky, which we estimate from Figure [7] as 4" 
(3 x 10 8 cm). Again, this is an oversimplification, but it sug- 
gests values of nu > 10 10 cm -3 within this target filament 
mass, comparable with those found for the much larger Tar- 
get 1. 

4.3. Comparison of methods 

Since the monochromatic method is more appropriate to 
fast-changing/-moving features, where we cannot co-register 
images taken at different times, it is of interest to gauge how 
well this method perfor ms wh en compared with our polychro- 
matic approach (Section [XT] ) for Target 1. 

Figure|9] shows correlation plots of monochromatic A/h (A) 
to polychromatic N^juted- ln the two shortest- wavelength 
channels, the signal to noise is rather poor, and the scatter 
in correlation is large. However, for the 171, 193 and 211 A 
channels, the scatter is less, and Nn(h) is shown to be con- 
sistently underestimated when compared with Abutted- In the 
same figure, we have colour-coded the points to show the fit- 
ted value of G; in general, the smaller the value of G, the 
greater the underestimate of the column density when only a 
single wavelength is used. This is perhaps unsurprising, since 
we set the filling factor arbitrarily to unity in calculating the 
monochromatic values. 

5. DISCUSSION 

Each of our approaches to determining A/h depends upon 
determining its value through the expression in Equation]?] 
We assume a fixed relati ve abundance for helium (Section |3f, 
but it should be noted that Gilbert et al. ( 2007]) and Kilper et al. 



( 2009 ) find variation of the He/H abundance ratio m promi 
nences, pointing to cross-field diffusion of neutral H and He 
as the most likely cause of stratification of this ratio. However, 
their finding is unlikely to apply since the dynamic timescale 
in the targets analysed here is m uch shorter than the diffu- 
sion timescale ( Gil bert et al.|2007| ). In eruptive events, [Kilper 
|et al.| ( [2009] ) also find that the ratio tends to homogenise in the 
hours leading up to eruption, and we suggest that the mixing 
of plasma induced by the eruption is likely to contribute fur- 
ther to this homegenisation. Work by Bempor ad et al.|(f2 009 ) 
uses the Hinode EUV Ima ging Spectrometer ( jKosugi et al. 



reduced field-of-view shown. However, we once again take 
profiles of Nn through a concentration of absorbing material, 



2007 ; Culhane et al. 2007 ) to show significant, spatially un 
resolved velocities (^ < 100 km s _1 ) at the location of an 
erupting prominence, which may support the picture of mix- 
ing. 

Since the absorbers of EUV radiation we consider here 
are ions and neutrals in a partially ionised plasma, we can 
only measure the total column density of neutral and par- 
tially ionised helium (thus, hydrogen), and estimate n^\ we 
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Figure 4. Left: Best- fit column density map for hydrogen, obtained by fitting measured absorption depth d (Equation|9}, to the function F (eqs.[lO|&[TT}. Right: 
as left, but showing the best-fit map of geometric factor G = fl^ /Iq. The field of view shown is identical to that shown in Figure[T] (See the electronic edition of 
the Journal for a color version of this figure.) 
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Figure 5. a) Total hydrogen column density, and inferred volumetric density, 
along the slices through Target 1 shown in Figure fflx b) Along the same 
slices, the value of the geometric factor, G is plottea! For both figures, the 
dashed line corresponds to the dashed slice in Figure^), while the solid line 
corresponds to the solid slice (See the electronic edition of the Journal for a 
color version of this figure.) 



can make no comment on the density of free electrons in this 
plasma beyond the trivial < n e < (hh/0.83), since the ratio 
n e /nu will depend upon the detailed parameters of the plasma 
and radiation field. Any fixing of the electron-to-hydrogen ra- 
tio to the usual coronal value of 0.83 is likely to lead to over- 
estimation of n e in the targets considered here. 

Similarly, we are unable to make estimates of the ionisa- 
tion degree of helium, since: 1) as we have shown, the cross- 
section for absorption of He I is very close to that of He II 
(Section [3]); and 2) we have no usef ul data longward of the 



He II ionisation edge (227 A; § [3JJ). Nor do we have data 
in both the helium continua and the H I Lyman continuum 
(504 < A < 912 A), which precludes measurements of He/H 
abundance variations. 

We also do not consider the information in the distribution 
of Nn about how thick Target 1 might be at each pixel - the 
concentration of mass in the lower parts is ignored, along with 
information that that may impart. The method can of course 
be inverted, assuming a constant value for nu to infer a spa- 
tially varying path length through the material. 

The approaches used in this article lack an advantage of the 
spectroscopic approaches favoured by previous authors. For 
example, we do not analyse emission from individual emis- 
sion lines, some with very short scale height, to make es- 
timates of the effect of lowered emission due to the cavity 
around a filament as Heinz eTet al.| ( |2003| ) do. However, in 
the case of an erupting filament, the cavity rapidly expands, 
so that both Iq and If are sampled within the cavity at the 
time of our measurements. Our use of emission from filters 
rather than at individual wavelengths is not a difficulty, how- 
ever, since the cross-sections for photo-ionisation depend only 
on wavelength. 



An earlier , spectroscopic analogue of the work by Kucera 
et al.| ( |1998]) applied to erupting prominences can be found in 
Penn (2000), who uses five emission lines between 555 and 
630 A to derive column densities of around 5 x 10 17 cm -2 . 
Although there are already likely to be significant differences 
between the amount of material studied by |Penn| and that stud- 
ied here, they also use lines in a regio n of the spectrum which 
has been shown by |Mein et alT1 < |2QQl| ) and |Gilbert et al.| ( |2011 ) 
to suffer from saturation of the photo-absorption effect. The 
much higher column densities we find may be due to the lack 
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Figure 6. a) Time average of AIA 193 A data between the times indicated above, representing the model of unattenuated intensity from the corona, Iq. b) Image 
from the same channel of AIA showing in-falling filament material that reaches what appears to be a magnetic null point, near position (560", -410"), at which 
the material is subsequently deflected east and west. In this image, the background radiation is attenuated in p laces by the filament material, and so represents 
lobs. c ) Column density of hydrogen, Nn, estimated with the monochromatic method described in Section [T2| The intersecting black and white lines at (550", 
-410") indicate the position of profiles of nn, shown in Figure [7] (See the electronic edition of the Journal lor a color version of this figure, where it is also 
available as an MPEG animation.) 
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Figure 7. Hydrogen column density, and inferred volumetric electron den- 
sity, along the slices through Target 2: the black/white line with la envelope 
corresponds to the black/white profile position shown in Figure^. (See the 
electronic edition of the Journal for a color version of this figure/) 

of these saturation effects at the shorter wavelengths studied 
here. 

5.1. Underestimation factors 

In Figures 15] and [6j we infer the volumetric electron density 
n e on the rignt-hanaaxis of the plotted profiles. In each case, 
we infer the density by estimating a path depth equivalent to 
the width of the feature observed. For the profiles in Target 1 , 
we estimate a path length of 14" (~ 10 9 cm); for the narrower 
feature in Target 2, we estimate a path length of 4" (~ 3 x 
10 8 cm). These depths are estimated from the FWHM of the 
corresponding A^h profiles. 

Studyi ng the He I D3 and H a lines in a number of promi- 
nences, Bommier et al. (1994 ) find mean volumetric electron 
densities of (2.1 ± 0.7) x 10 cm -3 . The values we infer 
from a monochromatic measurement of our targets are con- 
sistent with this mean, if a little on the low side. However, 
the monochromatic estimates are systematically lower than 
those derived from the multi- wavelength fit to d(X) by as 
much as an order of magnitude. In addition to this effect, 
there is the consideration that our estimates are again under- 
representations of the true column depth if a non- vanishing 
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Figure 8. Comparison of the best-fit H column depth in Target 1, derived 
polychromatially (top left), with that derived from individual wavelength 
opacities (remaining plots). Values are plotted logarithmically according to 
the colour scale in the centre. (See the electronic edition of the Journal for a 
color version of this figure.) 



fraction of the helium along the line of sight is fully ionised, a 
distinct possibility if heating tak es place in or around thi s ma- 
terial. In an erupting structure, Kucera & Landi (2008 ) find 
evidence of emission at log T e ~ 5.0. 

As already indicated above, the column depth can become 
so large that it effectively blocks all emission from behind. 
Gilb ert et"aL] ( |2011| ) recently encountered problems with the 
A 625 line, finding that the esimate of Nu at A 195 was much 
larger. The saturation of the effect at longer A will cause the 
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Figure 9. Correlation plots of the best-fit H column depth, on the horizontal 
axis, against that derived from monochromatic absorption. Colour-coding 
represents the corresponding geometric factor, G for each data point. With the 
exception of Nu(94) and A^h(131), the two lowest signal-to-noise channels, 
the monochromatic result is consistently an underestimate when compared to 
the fitted Nr. (See the electronic edition of the Journal for a color version of 
this figure.) 



depth to be underestimated: at some limiting value of T, no 
additional photons can be scattered/absorbed by the absorb- 
ing mass. Further evidence of this effect is shown by Heinzel 



et al.| ( |2008| ) as disagreement between estimates of the con 



tinuum opacity derived at Mg X A 625 with that inferred from 
THa. The former shows an underestimat e of column density 



compared with estimates from the latter. Gilb ert et al. (2011 
Ea ch paper attributes the problem of underestimated mass (as 
do Kucer a"et al.|1998| ) to the saturation effect described above. 

These effects would suggest that the ratio between Nh and 
the column depth of all absorbers deviates from the lower 
limit of Equation [4] However, there is an additional term 
which may lead to our underestimation of Nh- As stated in 
Section|3] we do not consider the effects of emission from the 
prominence material in the EUV passbands analysed. There 
could, in fact be emission in l ines over a wide range of temper- 
atures; O'Dw yer et al.| ( [20T0| ) have recently performed a de- 
tailed analysis of the emission lines (and corresponding tem- 
peratures) to which AIA is sensitive. Differential emi ssion 
measures construc ted from sections of a prominence by |Par-| 
|enti & Vial ( 2007) suggest that there is a considerable amount 
(though not a majority) of material at log 7^, > 4.8. Either of 
these cases would suggest that the effect on our results will 
be to make them a lower limit to the local filament mass, as 
already expected. Furthermore, in respon se to the work of 



Kucera et aL] ( |1998| ), |Engvold et aT] p001| ) point out that the 



former's results are susceptible to emission from hot lines in 
and around filaments before eruption. They are unable to es- 
tablish the contribution of these hot lines to the total observed 
emission, but filaments are certainly susc eptible to heating 



before and during eruption, as shown by |Kucera & Landi 
P008 1 " 
(e.g 



^008) , often with apparently coronal-tempera ture emission 
^.g., | Liewer et al.|[2009| |Sterling et al.||20lT] ) as is perhaps 



seen in Target 2 (Figure [6]). It is less clear wh ether there is 
heatin g of Target 1, the target in which we apply Kucer a et al. 
( 1998)'s method, as it returns to the lower solar atmosphere, 
but this possibility is certainly worth bearing in mind. In any 
case, the measurements made here would still serve as a lower 
limit. 

5.2. Advantages to this method 

Despite the disadvantages mentioned above, our ap- 
proaches have a number of advantages: 

1. They require no co-ordinated observing campaigns, 
only cotemporal or near-cotemporal images of an erupt- 
ing (or otherwise rapidly moving) filament in several 
EUV wavelengths. Full-sun EUV observations are car- 
ried out in a patrol like fashion all the time by AIA, so 
that the approaches can be applied whenever data are 
available. 

2. Because we use lines in the range where the helium 
cross-sections are at least as important as that of hydro- 
gen, the method is sensitive at higher temperatures than 
those which use the Lyman continuum only (504 < X < 
912 A). 

3. In exploiting the similarity of the He I and He II photo- 
ionisation cross- sections, we are able to measure the 
presence of all but bare helium, which allows us to es- 
timate the total hydrogen column density, not only that 
of neutral hydrogen, through an abundance argument 
(Equation]?]). 

4. The images used here are co-temporal or nearly co- 
temporal, allowing us to combine EUV image data in a 
way which has previously only been possible with spa- 
tially rastering or fixed- slit spectrometer data. 

5. When applied to erupting filaments in an expanding 
cavity, the complication of emissivity blocking is re- 
moved, as the cavity rapidly expands to encompass the 
observed area including To (but see below for a com- 
ment on this). 

6. Material from an erupted filament is likely to also be 
several scale heights into the corona, meaning that the 
majority of emission will be from behind the erupting 
material, as suggested by Figure]?] 

7. The method is, in principle, also applicable in a spatial- 
interpolative way. 

A qualification of point [5] is that emissivity blocking can 
be understood as low-temperature material replacing higher- 
temperature plasma that would otherwise contribute to the 
emission seen by AIA. In this sense, it is not limited to promi- 
nence cavities, and could in fact lead to some over-estimation 
of the absorption due to the filament, since the absorption 
depth is a measure of the decrease in intensity. However, since 
the filament material we have studied is likely to be rather 
high in the atmosphere, this consideration is unlikely to be 
i mportant here. 

|Heinzel et al.| ( |2008| ) draw an interesting (but probably in- 
correct) distinction between two types of observation previ- 
ously used to carry out this type of analysis: those typically 
only obtained at one wavelength by large-field of view instru- 
ments ("patrol-type"); and those that use a combination of as 
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many spectral lines as possible, termed "multi-wavelength" 
ob servations, largely pe rfor med by a super s et of t he authors 



Heinzel et 



al] (2008) and jKucera et al.] ( [T998] ). Perhaps 
Hack of co-temporal multi-wavelength data at 



because of the lack of co-temporal 
the time, their definition of "multi- wavelength" refers to spec- 
trometer observations, rather than images. We argue that the 
terminology is no longer valid since our method, which we 
term "polychromatic opacity imaging", uses the patrolling 
advantages of a full-time, full-Sun imager, in combination 
with the (near-)co-temporal images from multiple passbands 
in the photo-ionisation continuum to construct column density 
maps of large-scale features. Whereas contemporary EUV 
spectrometers are nearly always rastering, the examples we 
present have spatial simultaneity, and can therefore be applied 
to large-field observations. 

The results of fitted Nu shown in Figure [8] reveal a surpris- 
ing amount of coherence in their detail, ana this may reflec t 
agreement with a prior test carried out by Gilb ert et aLlp005| ). 
These authors compared the results of spatial- and temporal- 
interpolative approaches where possible, and found that the 
temporal-interpolative approach yields more accurate (and 
precise) results when compared with Thomson- scattering es- 
timates. 

Lastly, our approaches have been restricted to narrower 
fields-of-view in this article, but this is largely driven by the 
desire to test the techniques. In principle, the methods can be 
applied to full-disc images. SDO/AIA is well suited to this 
task, and the current rise toward solar maximum would sug- 
gest many suitable targets for the method. 

6. CONCLUSIONS 

We present a study which shows the powerful application 
of polychromatic opacity imaging to cotemporal and near- 
cotemporal data from SDO/AIA of a spectacular filament 
eruption. The single- and multi-pas sband approaches used 
in this study give consistent lower limits on the total hydro- 
gen column density (as opposed to just the neutral hydrogen 
column density) in the temperature regime below full ioni- 
sation of helium, and therefore are sensitive to a large range 
in prominence temperature. The values of column hydrogen 
density that we find are larger than in previous studies, but 
one of the targets chosen is rather exceptional and large col- 
umn masses are perhaps to be expected. 

Since we restrict our study to estimating the column density 
and inferring the volumetric density of hydrogen, this method 
removes the need for co-ordinated observing campaigns to 
study prominence mass through comparison of opacity at 
multiple wavelengths. Using AIA data allows us to perform 
such studies of erupting filaments wherever they are observed. 

Such observational constraints may be useful in the future 
development of more realistic models of erupting filament 
mass in coronal mass ejections. 
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